Cellular Ability to Sense Spatial Gradients in the Presence of Multiple Competitive 

Ligands 



Shu-Hao Lioi0 and Chia-Chu Cherjll 

Department of Physics, National Cheng Kung University, Taman, Taiwan 70101 

(Dated: December 26, 2011) 

Many eukaryotic and prokaryotic cells can exhibit remarkable sensing ability under small gradient 
of chemical compound. In this study, we approach this phenomenon by considering the contribution 
of multiple ligands to the chemical kinetics within Michaelis-Menten model. This work was inspired 
by the recent theoretical findings from Bo Hu et al. [Phys. Rev. Lett. 105, 048104 (2010)], our 
treatment with practical binding energies and chemical potential provides the results which are 
consistent with experimental observations. 
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I. INTRODUCTION 

Cellular sensing ability is a general but critical biolog- 
ical function, it plays important roles in cancer sensing, 
wound healing, embryogenesis and neuronal development 
This remarkable ability allows cell to gain necessary 
energy and nutrition and obtain information from other 
cells Ingeneral, the size of the cell is only a few mi- 
crometers [S] but it can discriminate the correct direction 
of the gradient when tiny variance of chemical concen- 
tration. The minimum fluctuation about detecting the 
concentration can be seen as the accuracy of the sensing 
ability, which was addressed by Berg and Purcell ^] and 
then modified by R. Endres et al. who considered the un- 
occupied time intervals of the ligands and receptors Q. 
Moreover, other researchers reported the physical limits 
of spatial sensing ability 0, and then demonstrated the 
impossibility to increase the ellipse cell's sensing ability 
by enlarging the cell's body Q- 

According to the report from Bo Hu et al. [6], the 
accuracy to sense the ligand gradient direction will be 
increased dramatically by larger cell. Their results are 
based on the basic assumption that the thermodynamics 
equilibrium and the chemical kinetics are equivalent. The 
thermodynamic results are obtained in the frame work of 
canonical ensemble. The gradient sensing ability is calcu- 
lated from the partition function which is completely de- 
termined by the energy of the Hamiltonian. The results 
of are very interesting and intriguing in the sense that 
the sensing ability can be established in such a simple 
model. In this work we modify the model of Q such that 
the correlation between energy and concentration can be 
relaxed. This is done by considering the dynamics of lig- 
ands. The chemical equilibrium of ligand and receptor 
provides a way to break the energy-concentration con- 
straint by introducing chemical potentials for addressing 
the problem of concentrations. In this work we have re- 
analyzed the problem of sensing ability with the grand 
canonical ensembles with the contribution of ligands in- 
cluded. Furthermore it is also interesting to address this 
sensing problem in a more general environment in dif- 
ferent ligands. In our model, the cell can distinguish 
distinct ligands with different chemical dissociation con- 



stants. This result can only be achieved by considering 
the binding energies of both ligands. Fortunately, our 
approach can also be extended to the analysis of this 
general environment. 

In this report, the assumptions and results of Q are 
discussed in section |lTl Our approach will be presented 
in section IIII[ and we will have a brief discussion on 
the experiment results related to our model. Section ITVl 
provides the analysis of multiple ligands system within 
Michaelis-Menten model. The conclusion is given in sec- 
tion IVl 



II. APPLYING ISING MODEL IN SENSING 
PROBLEM 



To deal with the cellular sensing problem, Bo Hu et 
al. treated the receptors as an Ising spin chain. From 
this model, they could derive the asymptotic variances 
as functions of the gradient steepness p and direction of 
gradient (j) which are parameters related to the concen- 
tration. According to the Cramer-Rao inequality, these 
variances determine the lowest uncertainties of cell sens- 
ing ability [8]. Here we briefly review their calculation 
procedures. 

In this model, the cell with diameter L and N receptors 
were immersed in the concentration environment which 
contains identical ligands. In their work all results were 
calculated with N = 80000, which is close to the prac- 
tical situation Q. The local concentration of ligands at 
nth receptor is Cn = Co exp[| cos((^„ — (j))], where Co 
is the background concentration, p is the steepness of 
gradient [p = ^|vC|), (p„ — 2mT/N denoted the loca- 
tion of Tith receptor and direction of gradient is (/). In 
this approach, ignoring the dynamics of ligands, the sys- 
tem is completely described by receptors which have only 
binding state {Sn = +1) with energy — e„ and unbind- 
ing state {Sn = —1) with energy e„, where e„ is given 
in the unit of thermal energy ksT . Due to the Boltz- 
mann distribution, the binding probability of the nth re- 
ceptor is Pon = jF7q:^^-r;r- Using simple receptor- ligand 
kinetics, the binding probability of the nth receptor is 



where Kd is the dissociation constant. By 
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FIG. 1. (Color online). The Cov[z\, Z2] of both models under 
Kd = Co and (fi = n/5 arbitrarily. The dashed line repre- 
sents Cov[zi, Z2]/a-^ versus p in [0], and solid line represents 
Cov[zi,z2] with p in our model which will be de- 

y<A(zi)2><A(z2)"> 

scribed in next section. 



FIG. 2. (Color online). The diagrams of our model: The 
cell with identical receptors is located in concentration pool 
with particular direction of gradient. Each receptor can only 
sense the ligands inside the sensing volume v. Three different 
energy levels with binding, unbinding and ligand states are 
all independent of location. 



assuming these probabilities are identical, the free energy 
is 



1 , Co ^ P , 
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By defining three statistical quantities (20,21,2:2) — 
iJ2n^n,J2n ■^Sn COS ipn,J2n ^ sin ) and transforma- 
tion {ao,ai,a2) = (^ In ■^,pcos(^,psin(/()), the Hamil- 
tonian is given as iJ^rjS'n} — — ao-^o ~ ^i"i+^2'a2 ^^^^^ 
by computing the logarithm partition function In Qjv = 
7Vln(2coshao) + g^^^.^'^^ + OCp"), they obtained the 
expectation values and fluctuations of zi and Z2' 
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Under the assumption Cov[zi, 22] = @ the joint prob- 
ability density of 21 and 22 is the Gaussian function, the 
asymptotic variances Up and cr^ can be obtained as: 
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These asymptotic variances are the minimum fluctua- 
tions with unbiased estimation of p and (f> vifhich are re- 
lated to sensing ability. Therefore, the fluctuation ct^ 
drops as gradient increase, and the sensing ability will 
be increased with enlarged cell volume since p on L and 
N (xL^ with 0<S <2 J^. 

In passing we would like to show the detailed analysis 
about the covariance Cov[zi, Z2] — for justifying the 
joint probability density as Gaussian distribution. The 
Cov[zi, 22] is given by: 
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and the numerical results is plotted in FIG. [I] in which 
one can see that it is very small and can be approxi- 
mated to zero, such that the probability density can be 
described by Gaussian distribution. Even if our model 
will use this approximation for further analyzing and 
Cov[zi, 22] in our model is also shown in FIG. [T] 



III. ADDING THE LIGANDS INFORMATION 

To calculate the physical limit of multi-ligands system, 
the ligand's concentration should be considered in the 
partition function. The grand canonical ensemble is then 
appropriate to construct our model. Here we adapted the 
notations of |6( where the cell with size L has N receptors, 
and the spatial information of ligands are given by the 
concentration C„ , gradient steepness p and the direction 
of the gradient 0. Each receptor should only sense the 
ligands inside an identical sensing volume independently. 
In other words, the nth receptor's sensing volume v is 
completely separate with others'. Let L„ = vCn denotes 
the number of ligand which can be sensed only by the 
nth receptor. The approximative expression of sensing 
volume is w = Airr'^d/N, where r — L/2 is the radius of 
cell and d is the size of ligand. However, to be consis- 
tent in this model, three different energy levels were set 
up to describe different states: the unbinding energy £„, 
binding energy Sb and ligand energy and the corre- 
sponding chemical potentials are denoted by fib and 
fii respectively. We notice that instead of being position 
dependent as resulted in Q, these energy levels are po- 
sition independent as they should be according to basic 
quantum principle. The Hamiltonian of this system is 

71=1 



(7) 



To simplify our discussion, we denote the nth recep- 
tor's grand canonical partition function as the binding 
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part Zbn and unbinding part z^n, which are given by: The factorial factors appear here since ligands are all 

identical. The total grand canonical partition function 

^, _ r,-/3f(et,-M6)+(£i-Mi)(-L„-i)] (o„\ of the whole cell becomes 

^""^ (i„-l)!' ' ^^^> 
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The binding probability due to the Boltzmann distri- 
bution is Pbn = -^f^^ — • The binding probability of 
chemical equilibrium at the nth receptor is given by 
Pen = JT^^K'- imposing Pfc„ = Pen we obtain the 
relation: 

a=^^ ^_g-/3[-(£i,-A'6) + (e,.-A'u) + (£i-A<i)], (10) 

Co wCo 

Eq. ([T0|) shows that the chemical dissociation constant 
Kd under this assumption should depend on energy lev- 
els, chemical potentials and the sensing volume. More- 
over, the binding probability would depend on concentra- 
tion and position, but the energy levels are independent 
of position. By using Eq. pH)) all the dependence on e 
and [i are replaced by dissociation constant Kd and local 



concentration Cq. In particular, the sensing ability will 
be determined by a = Kd/C^ and more details will be 
discussed as the following. 

With the statistical quantities (zo,zi,Z2) = 
{Y.,i^n,Y.nl'^ri cos (p,i,Y.nl^ri sin (fn) and the trans- 
formation (q!i,Q!2) = (p cos sine/)), the analysis 
of discrimination can be proceeded with the expec- 
tation values 'zi and Z2. These expectation values 
are calculated by using direct summation method. 

Defining C„ = e*'^", then zT = Rei^n Cn [lll+lZ] ^ ^^^^ 
= Im^^n Cn (^|^"+z""j ]- Under the small p assump- 
tion (which is also true for experimental environment), 
we expand the summation to second order of p, and 
treat the summation as an integral over [0, 27r], hence 
the integrals can be computed as 
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and both results of zi = Re[n] and Z2 = Im[n] are: The fluctuations of zi and Z2 are =< A(zi)^ >= 

i?e[cr^] and cr| =< A(z2)^ Re[a^] where 
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We should mention that all the results are real, al- 
though complex numbers appear from complex integral 
calculations to simplify the notation. It is interesting to 
note that even though the grand canonical ensemble de- 
pends on energy level e and chemical potential /i, in this 
model the expectation values and fluctuations only de- 
pend on p, (j) and a, where a implicitly depends on e and 
/i as given in Eq. (jlOp . The detail dynamics of the system 
such as the energy levels and the chemical potentials de- 
termine the dissociation constant Kd- However, the sens- 
ing ability relies on the detail dynamics only through the 
specification of a. In Q they used the procedures given in 
[8| to estimate the joint probability as Gaussian function. 
To apply this assumption one should check the value of 
Cov[zi, Z2] explicitly. By using Eq.® the correlation 
which was divided by fluctuation is obtained in FIG. [TJ 
where one can see clearly that Cov[zi, Z2] ~ as p being 
small. Under this assumption, according to [8] one can 
assume the joint probability density as Gaussian func- 



tion, i.e., /(zi,2|zi,2) ~ 2^^exp[ ^^.^ 

However, for large p, it seems that a more detailed analy- 
sis is called for. Under the Maximum Likelihood Estima- 
tor (MLE) theorem, the estimators of p and (j), denoted 
by (Tp and cr^ respectively, can be obtained near the ex- 
pectation value without any bias Q. The estimators ( 
also asymptotic variances) related to the sensing ability 
can be derived by the Cramer-Rao lower bound: 
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FIG. 3. (Color online). The direct summation results of ap 
(black) and ct^ (orange) versus p with a = 1.0 or Ko — Co- 
The dash line represents the results in [^. The triangle rep- 
resents the (T^ under a = 0.01, which means the local concen- 
tration is much larger than the dissociation constant Kd- 



decrease when the steepness increases under the condi- 
tion Kd = Co- We also plotted their results [g] in FIG. [31 
one can see that the result of two models are quiet close 
with similar characteristics. 



It is interesting to see our results are related to experi- 
ment observations. It is well known that the sensing abil- 
ity or sensing accuracy strongly depends on the steepness 
of gradient but weakly correlates to background concen- 
tration 10]. In this work we have analyzed this aspect 
and the results are also plotted in FIG.[31 cr^ approaches 
to zero when p is large, which indicates good sensing 
ability in this range. When p is larger than 0.1, it can 
be seen that cr^ of different concentrations (a = 0.01 
and a = 1000) are more or less the same, which means 
that the local concentration will weakly depend on sens- 
ing ability when p > 0.2. However, when p is less than 
0.1, the effect of concentration might be large for sensing 
ability. 




Following figures are the computational results with- 
out any approximation. FIG. [3] shows that ct^ which 
relates to the ability of sensing will decrease by increas- 
ing p. Since the steepness of gradient p is proportional 
to the size of cell L, the sensing ability would increase 
dramatically when volume of the cell expands. The re- 
sult shows similar conclusion as @ without £„ being C„ 
and position dependent. In FIG. [31 one can see that 
almost remains the same for different p, however cr? will 



Moreover, it is known that cells can show remarkable 
sensing ability in particular steepness and concentration 
range. For instance, Dictyostelium cells will move to- 
ward cyclic adenosine 3',5'-monophosphate (cAMP) to 
function as a chemoattractant [111 ]. In general, under 
p = 2% and Kd ^ lOOnM, the cell can exhibit the sens- 
ing ability when cAMP is in a range between lOpM to 
lO/xM dl], or a = 0.01 - 1000 in our model. FIG. 
[4] shows the two asymptotic variances ap and cr^ within 
this range. It is obvious that the fluctuation is quite large 
at a = 1000 (which means small background concentra- 
tion or large disassociate constant), but the cr^ is below 
1.0 when p > 0.32. In other words, the sensing ability is 
exhibited under small concentration and small steepness 
of gradient. 
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FIG. 4. (Color online). The results of cTp (black) and tj^ 
(light orange) versus p in different background environment. 
The dash line represents a — 0.01 and the solid line shows the 
result of o = 1000. The result is plotted by direct summation 
without any approximation. 



IV. MULTIPLE LIGANDS WITH 
COMPETITIVELY BINDING IN CELLULAR 
SENSING ABILITY 

In practical biological systems, many receptors can 
bind with different ligands on the same site to conduct 
many important functions [T3]. For example, the ions 
H^, K"*" and Mg^^ can bind with eukaryotic cell's Ca'^'^- 
binding sites of calmodulin [l^l , which are related to the 

I 



intracellular movement, metabolism and apoptosis [15|. 
Moreover, different ligands might present dissimilar ef- 
fect after binding. For instant, platelet-derived growth 
factor (PDGF), which is excreted by platelet a-granules 
during the injury, can strongly attract monocytes and 
neutrophils [16[. However, the protamine sulfate which 
competitively binds to the surface of monocyte and neu- 
trophil will block the chemotaxis |l6t] and shows distinct 
role in chemotaxis. Therefore it is interesting to analyze 
the sensing ability in the multiple ligands environment. 

In this section we consider the concentration environ- 
ment with two different ligands denoted by Ligand 1 and 
Ligand 2, which can bind on the same site of the receptor 
whenever it is unoccupied by any ligand. The receptors 
are still be treated as Ising spin chain. Since the ef- 
fects of binding for these ligands are in general distinct, 
two different operators Sin and S2n are needed to de- 
scribe the states of receptor for Ligand 1 and Ligand 2 
respectively. The eigenvalues of both operators are equiv- 
alent, such as binding state (-1-1) and unbinding state 
(—1). The parameters of two ligands are listed in TA- 
BLEU where the concentrations of both ligands are C„ = 
Coexp[^cos((/?„-(/)i)] and Ai = ZJq exp[^ cos(i^„-02)] 
at the nth receptor, and the energy level and chemical po- 
tential for unbinding receptor are e„ and fiu respectively. 
All the energy levels and chemical potentials of the sys- 
tem are independent of position and concentration. The 
Hamiltonian of the system is expressed as: 
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V fC CI Y^r / ('5'ln + l)('S'2n - 1) X , / ('S'2ri + 1 ) (S'ln ~ 1 ) s, , / (-Sin - 1) (5'2n " 1) x 
nN{t3in,i32n\ =2^[^l 7 ^ ( 7 ^ ( 7 i^u 



r ^ {Sln + l){S2u-l) ,^ _U r n ^ (^2n + l)(^ln-l) ,„ . 

+ (w6„ H )£i + [vDn H )e2\ 



(16) 



TABLE I. The parameters for two ligands at the nth receptor. 



Type 



Ligand 1 Ligand 2 



Background Concentration 


Co 


Do 


Steepness of Gradient 


Pi 


P2 


Direction of Gradient 


01 




Energy Level of Unbinding Ligand 


£i 


£2 


Energy Level of Binding Ligand 


e\ 




Chemical Potential of Unbinding Ligand 


Ml 


^J.2 


Chemical Potential of Binding Ligand 


Ml 





The grand canonical partition function at the nth site 
Zn can be separated as parts of unbinding (z„„), binding 
with Ligand 1 {zin) and Ligand 2 {z2n)- 



1 



{vCnY-ivD^y. 



,-/3[(e„-/.„)-H(ei 



-/il)«C„ + (£2-M2)t'-D„] 

1 

(17a) 



1 



zi 



(vCn - ly.ivDny. 

g-/3[(£'i-A<'i) + (ei-Mi)(t'C„-l) + (e2-M2)i'£'„]^ (17b) 



Z2: 



{vCny.ivDn - 1)! 

g-/3[(£2-^2) + (6i-Aii)uC„ + (e2-A'2)(f£',i-l)]^ (17c) 

The total partition function S for the cell therefore be- 
comes 



AT AT 
5 = J]^ Z„ = J]^ {Zun + Zin + Z2n)- 
n=l 



(18) 



n=l 



To identify the probabilities of binding and unbinding 
cases in chemical kinetics, we apply the Michaelis-Menten 
Model which is widely used in non-allosteric enzymes [17 1 
where the receptors are treated as isolated enzymes in 
this model. The parameters of Michaelis-Menten Model 
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TABLE II. The parameters of Michaelis-Menten Model at the 
nth receptor. Abbreviation: IC, initial concentration; FC, fi- 
nal concentration after equilibrium; CDC, chemical dissocia- 
tion constant. 



Type 



Parameter 



IC of Enzyme 

FC of Unbinding Enzyme 

FC of Enzyme Bind with Ligand 1 

FC of Enzyme Bind with Ligand 2 

CDC for Ligand 1 and Enzyme 

CDC for Ligand 2 and Enzyme 



en 

Cnl 
C„2 

K2 = fcj/fe 



are listed in Table |TT1 The main equations for competi- 
tive constrain and chemical equilibrium for Ligand 1 and 
Ligand 2 can be written as: 



Cnl + C„2, 



(19a) 



dt 



klCnCn - hCnl = 0, 



(19b) 



dc„ 



dt 



fc2e„i?„ - k2Cn2 = 0. 



(19c) 



At equilibrium, one can solve for the final concentra- 
tions Cnl, c„2 and the unbinding receptors. The ratios of 
these variables and the original concentration of receptor 
can be seen as the probability of binding or unbinding 
state at nth site. Therefore, the probabilities of binding 
with Ligand 1 (Pin), Ligand 2 (P2n) or without binding 
{Pun) are: 



e„o l + §7 + §7 ^« 1 + l,C'„e-'3[('^i-A'i)-(el-/'l)-(e"-M^.)] + „£)^je"'3[('^2-A'2)-(£2-A'2)-(£„-M„)] ' 



Ki 



Zlr 



y(7^g-/3[(ei-M'i)-(ei-A'i)-(eu-Mu)] 



p ^ = 

^" e„o 1 + %^ + ^ Zn l + uC„e-'5[(^i-^i)-(^i-^i)-(^"-*'")l +t;D„e-^[(^2-''2)-('^2-A^2)-(e.-A'J] 



(20a) 
(20b) 



P2« = 



Cn2 

SnO 



K2 



_ Z2n 
C„ , Z5„" = ~y~ 



u£'„e~''''^2-^'2)-(^2-A'2)-(eu-A't.)] 



t;C„e~'^[(^'i~'^'i)~(^i~''i)~(^"~'^")l -I- tiD„e~^[(^2-M2)-(e2-M2)-(eii-/Jii)l ' 
I 



(20c) 



Similar to the case of single ligand system, we introduce 
h\ and h2- 



Co 



h - ^2 

02 — 

Do 



1 



vDo 



-/3[-(e;-Ai'i) + (ei-pi) + (e„-p„)] 



- /3 [- (^2 - /^2 ) + ('=2 - /i2 ) -I- (Eu - /Ju ) ] 



(21a) 



(21b) 



It is necessary to define the statistical pa- 
rameters for Ligands, which are {xi,X2) — 



(I]n ^'S'inC0S(y9„,X;„ ^-SinSin-^n), and (1-3, 2^4) = 
(J2nh^'2riCOSipn,J2n^^'2riSiii(pn) for Ligaud 1 and 2 
respectively. Furthermore, the expectation value for 
each parameter is = < Xi >, and it's fiuctuation 
can be calculated by erf = < (xi — Xi)^ >. Assuming 
these parameters are independent, their probability 
distributions can be described by Gaussian function 



nti2¥^exp[ 



- f ], then the Fisher informa- 



tion matrix F then can be obtained as 
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< 
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a" Ing 

dxidx2 
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a^ Ing 

dxidxs 
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62:2^ 
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a" Ing 

dx20x3 
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a" Ing 

aa;2aic4 
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d'xi dx3 
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a" Ing 

dx2dx3 
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a" Ing 
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a" Ing 

dxsdxi 


> 
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a" Ing 
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a" Ing 

aa;2aa:4 


> 


< 


a" Ing 


> 


< 


a" In g 

aa;42 


> 


/ 



(22) 



I 

The asymptotic variances of 0^ (fc = 1,2 for two lig- ands respectively) can be obtained by the following in- 
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FIG. 5. (Color online). The results of cr^^ and cr^^ versus p 
for bi — 0.01 with different 62 with <j!>i = and <^2 = tt/IO. 
We just plot the range 0.02 ~ 0.2 due to the lack of space. 
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FIG. 6. (Color online). The results of a'^^ and a^^ versus p 
for 61 = 100 and various &2 with <j!>i = and <j>2 = tt/IO. 



equality 



'Rfe > 0, 



(23) 



where is the transformation vector of F is the 
inverse of information matrix. is the minimum fluc- 



tuation of ( 



The definition of Ri, is: 



R 



dxi 



dx2 0X3 &XA 



(24) 



By expanding the C„ and I?„ to second order of p, we 
obtain 6i and 62 as: 



arctan 



{2 + l-)x2-{l-i)x, 



bi 

(2 + ^)5^1 



(1 
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FIG. 7. The results of cr^^ and cr^^ versus p for fei — 100 and 
62 = 400, 4000 with <f>i ^ and 4>2 = vr/lO. 
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The minimum of Cfe is obtained by inequality [Eq. (j23p ]. 
By taking P2 = Pi = P, FIG. [5] and [5] are the results of 
fixed bi with different 62 by numerical differentiating Eq. 
(|25p with statistical parameters a;^. 

One would expect that small asymptotic variance 
presents better sensing ability of the ligand. Due to the 
lack of experimental evidence for cr^, we set a criteria 
< 1 to be the condition for sensing ability. In FIG. 
[3 with bi = 0.01, 62 = 0.4 and p = 0.02, cr^^ ~ 300 
and cr5 — 1000, it seems that the cell can not have a 
good sensing to any ligand. As p increases, for exam- 
ple p = 0.4, = 0.81864 and ct^^ = 2.50598, which 
means the cell has small asymptotic variances for better 
resolution of the ligands. One may infer that the cell 
might recognize a more distinct direction and move to- 
ward Ligand 1. Meanwhile, the difference of the asymp- 
totic variances are not so eminent when these bj (j = 1,2) 
are almost equal, this can be seen from FIG. [S] that for 
the case 62 = 0.02 and p = 0.61, ct^^ = 0.98659 and 
cr? = 1.01726, and the cell could not determine a prefer 



direction for motion, even though cr^ satisfy our criteria 
but the cell might be undecidable due to cr^^ ~ a^^ . 

It is interesting to see whether the above observation 
depends on the magnitude of bj. The results from fixed 
bi = 100 with different 62 are plotted in FIG. El For 
the case of 62 = 400, we have cr^^ = 14.99085 and 
= 187.88 at p 0.04, and cr^^ = 0.95067 and 
cr^^ = 11.91643 for p = 0.16. It seems that better sensing 
resolution exhibits with larger gradient. On the other 
hand, for the case 61 = 100, 62 = 120 and p = 0.2, 
cr^^ — 0.812736 and cr^^ = 1.1280, confusing information 
also exists in this situation. Therefore, FIG. |6] indicates 
the occurrence of good resolution with large difference of 
bj is a general feature in this theoretical model. How- 
ever, further experimental observation will be needed to 
clarify this assertion. 

Since large bj corresponds to unbinding, it is important 
to see the effect of unbinding on FIG. El For large 62, such 
as 62 — 4000, the results show the cell has no information 
for Ligand 2 through binding, therefore one would expect 
a large cr^^ and our results in FIG. [71 also shows such 
tendency. Furthermore, when 62 decrease to 400, the cell 
can have more binding receptors , therefore cr? decreases. 
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Above discussion indicated better sensing resolution 
exists when both conditions < 1 and large difference 
of (T^^ and (T^^ are satisfied. In addition, the asymptotic 
variance decreases when the receptor tends to bind with 
ligand and cell could receive more information from lig- 
ands. 

We have provided the extension to deal with two lig- 
ands system by using Michaelis-Menten model. Under 
the equilibrium situation, minimum fluctuations of sens- 
ing ability can be obtained. It is noted that such system 
can be analyzed by including the ligand concentrations 
inside the partition function. 

V. CONCLUSION 

In this work we have modified the mechanism of sens- 
ing ability by including the dynamics of the ligand. In our 
approach we are able to avoid having energies correlate 
with the concentration of ligand. This was accomplished 
by setting up the system with different energy levels and 
chemical potentials, and use grand canonical partition 
function to address the sensing ability. It is interesting 



that this model can still exhibit remarkable sensing abil- 
ity. Moreover, in our approach the energy levels are free 
parameters which can be used for any cellular complex, 
therefore our model has predicting power for other phys- 
ical quantities such that further experimental results can 
be used to justify the correctness of this kind of model. 

We have also studied in more complicate environment 
with multiple ligands under the competitively binding 
situation. The results indicate that cells can distinguish 
different ligands under the small gradient with different 
chemical dissociation constants. This extension gives the 
possible way to predict the reactions of the cell under the 
practical biological environment. 
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